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Abstract. - In statistical physics, the conservation of particle number results in the equalization 
of the chemical potential throughout a system at equilibrium. In contrast, the homogeneity of 
utility in socio-economic models is usually thought to rely on the competition between individuals, 
leading to Nash equilibrium. We show that both views can be reconciled by introducing a notion 
of chemical potential in a wide class of socio-economic models, and by relating it in a direct way to 
the equilibrium value of the utility. This approach also allows the dependence of utility across the 
system to be determined when agents take decisions in a probabilistic way. Numerical simulations 
of a urban economic model also suggest that our result is valid beyond the initially considered 
class of solvable models. 
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Qh Socio-economic sciences and statistical physics arc both 
interested in the evolution of systems characterized by a 
£Sj large number of interacting entities. These entities can 
for instance be economic or social agents in social sciences 
*^ j atoms or molecules in statistical physics [1HS] ■ The 
question of the emergence of macroscopic patterns from 
^{■^ the interactions of a large number of microscopic agents is 

'studied by both fields of science. In statistical physics, a 

^— ^ quantitative framework has been developed over the last 
f— ^ century, allowing the equilibrium behaviour of large as- 
i— H semblics of atoms or molecules to be handled precisely [B] . 
^ In socio-economic models, the preferences of individuals 
• ^ are usually characterized by a utility function, which de- 
scribes their welfare with respect to their current situation 
or environment. Each individual or agent wants to max- 
imize his own welfare. Decisions (e.g., moving to a more 
convenient place) are thus taken in a purely selfish way, 
while in physics the motion of particles is governed by the 
variation of the total energy. Recently, a global function 
linking individual decisions to the variation of a global 
quantity has been introduced to describe some classes of 
socio-economic models [HI!]- This approach then allows 
such models to be described with statistical physics tools. 
Importantly, the equilibrium state can then be calculated 
by maximizing a state function (akin to a free energy) in- 
stead of having to solve a complicated Nash equilibrium 
of strategically interacting agents. 



The question we investigate in this letter is whether this 
physical description of socio-economic models can be ex- 
tended to other basic concepts of statistical physics, such 
as the equalization of thermodynamic parameters like tem- 
perature or chemical potential. The equalization of these 
quantities throughout the system precisely results from 
the conservation of the conjugated extensive quantities, 
namely the energy or the number of particles. Although 
there is no notion of energy in socio-economic models, the 
dynamics indeed conserves the number of agents. A natu- 
ral question is thus to know whether a chemical potential 
can be defined in such models, and what would be its re- 
lation to standard socio-economic concepts. This question 
is further motivated by the following remark. In spatial 
socio-economic models, the individual dynamics leads to 
a Nash equilibrium, where no agent has an incentive to 
move. If all agents are of the same type, the Nash equilib- 
rium results in a spatially uniform utility, even if the en- 
vironment is spatially inhomogeneous like in cities, where 
the center plays a specific role. This uniformity property is 
also expected from the chemical potential (if such a quan- 
tity can be defined), suggesting a possible relation between 
these two notions. 

Here, we investigate this issue in the framework of a 
generic class of exactly solvable models involving a popu- 
lation of locally interacting agents. We define in a precise 
way a chemical potential for this class of models, and pro- 
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vide a direct link between the chemical potential and the 
socio-economic utility. Two explicit examples from the 
field of urban economics are also presented. 

Model and dynamics. — This work deals with socio- 
economic models characterized by a large number of in- 
teracting agents, residing on a set of sites, labeled by an 
index q = 1, . . . , Q. Agents are able to move from one 
site to another in order to increase their utility. In ad- 
dition, agents belong to to different groups, according for 
instance to their income, or to their cultural preferences. 
The variables used to describe the system are the numbers 
n q i of agents of each group i — 1, . . . , m at each node q. 
The configuration of the system is described by the set 
x = {n q i}. We assume that agents cannot change group, 
so that for all i, the total number Ni = J2 q n qi °f agents 
of group i is fixed. The satisfaction of agents of type i on 
site q is characterized by a utility U q i(n q \, n qm ) that 
depends only on the numbers of agents of each group on 
the same site q. 

The model is defined with a continuous time dynamics 
following a logit (or Glauber) rule, which is commonly 
used in social sciences and in particular economic works 
[5] . If transitions between sites q and q 1 are allowed, agents 
move from q to q' with a probability per unit time 



W = 



1 + e -AU/T> 



(1) 



where AC/ = U qH — U q i is the variation of the agent's own 



utility, with 

u qi 
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(2) 
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The parameter T plays the role of an effective tempera- 
ture, introducing noise in the decision process to take into 
account other factors influencing choices [5], and vq is a 
characteristic transition frequency. 

In order to obtain analytical results, we assume that the 
utility function is such that the change of individual utility 
experienced by an agent during a move can be expressed 
as the variation of a function of the global configuration 
x = {n q i} [7]. More precisely, we assume that there exists 
a function L(x) such that for each agent in group i, moving 
from node q to node q' , 



U a 



L(y) - L(x) 



(4) 



where y = (n u , . . . ,n qi -l, . . . , n q >i + l, . . . ,n Qm ) and x = 
(nil, ■ • ■ , TiQm) are the configurations of the system after 
and before the move respectively. Such a function L(x) 
thus provides a link between the individual behaviour of 
agents and the evolution of the whole system. In physical 
terms, it can be thought of as an effective energy. The 
relevance of this assumption (which bears some similarities 
with potential games [TU|) for the general class of systems 
considered above will be discussed at the end of the paper. 

The stationary probability distribution V s {{n q i}) = 
V s (x) is obtained by solving the master equation governing 



the dynamics of the system [TT]. If Eq. Q holds, detailed 
balance is satisfied [HJIH], and we obtain the following 
expression for the distribution V s (x): 



V s {x) = 



1 



,L{x)/T 



z s n 



q,i n qi 



(5) 



where Z s is a normalization constant. The product of 
Kronecker S functions accounts for the conservation of the 
total number of agents in each group. The different fac- 
tors appearing in Eq. ([5]) can be given a simple interpre- 
tation. The exponential factor directly comes from the 
detailed balance associated to the logit rule Eq. ([T]), while 
the product of factorials appearing at the denominator in 
Eq. ([5]) results from the coarse-graining of configurations. 
Namely, given the numbers of agents there are for 

each group Ni\/J\ n q i\ ways to arrange the agents of the 
group. As the numbers Ni arc fixed, Nil can be reabsorbed 
into the normalization constant. 

Defining a density of agents p q i = n q i/H, where H 3> 1 
is a characteristic number (for instance a maximal number 
of agents on a site) , the utility U q i then becomes a func- 
tion u q i(p q i, p qm )- We further assume that the function 
L{x) can be written in the large deviation form |13j 



L(x) = HL({p qi }). 



(6) 



To determine L, we combine Eqs. (|4]) and ([6]), and expand 
L to leading order in 1/H, yielding 



dL dL 



OPq'i OPqi 

By identification, we get for all q 
dL 



dpq 



Hqii^Pql-) •••■> Pqrn)' 



(7) 



(8) 



As the r.h.s. of Eq. © only depends on densities of agents 
on node q, L necessarily takes the form 



and one has 



L{{Pqi}) = ~Y^l q {p q \, pqr, 



dig 

dp 



(9) 



(10) 



If there is a single group (to = 1), l q {p q ) is simply obtained 
by integrating u q {p q ). In contrast, if m > 1, l q (and thus 
L) only exists if the following condition, resulting from the 
equality of cross-derivatives of l q , is satisfied: 



du qi _ du qj 



i ^ 3- 



(11) 



Opqj dp qi 

If this condition holds, the stationary distribution reads, 
after an expansion of the factorials using Stirling's for- 
mula, 

v{{ Pql }) = |n eH/9(pai '-' P9m)/T n<E^ - qp. 

q i q 

(12) 
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where f q is given by 

fq(pql, —,Pqm) = lq(pql, —, Pqm) + Ts(p ql , p qm ), (13) 

with 

s{Pql, Pqm) = -^Pqi^Pqi- (14) 

i 

In analogy to physical systems, f q (p q i, p qm ) can be in- 
terpreted as a local free energy (up to a change of sign), 
and the term s(p q i, p qm ), which is multiplied by the 
'temperature' T, may be seen as an entropic contribution 
associated to the node q. 

Utility and chemical potential. — We now turn to 
the main result of this letter. The configurations {p* qi ] 
which maximize F = ^2 f q under the constraints of fixed 
global density J2 q p q i — QPi are the most probable (or 
equilibrium) configurations. Finding the equilibrium den- 
sities of agents is then a constrained maximization prob- 
lem. Let us introduce a Lagrangian 

£({P«»}> {*i}) = ^2fq(Pql:->Pgm) ( i5 ) 
Q 

~ Yl Xi ( Yl ph - QPi) ' 

i q 

where the parameters Xi are Lagrange multipliers associ- 
ated to the conservation of the number of agents in each 
group. In physical terms, A^ corresponds to the chemi- 
cal potential of the agents of group i. The equilibrium 
densities {p qi } are then determined from the conditions 
dC/dp q i = for all (q,i), yielding 

ds 

Uqi(p ql ,-,P qm ) + T- {p q i,-,P* qm ) = K, (16) 

QPqi 

which is the main result of this letter. Equation (|I6|) thus 
provides an answer to the question raised at the beginning 
of this letter: there is indeed a direct relationship between 
the socio-economic utility and the chemical potential de- 
fined, in analogy to equilibrium physical systems, from 
the conservation of the number of particles. At zero tem- 
perature, both quantities can be identified. This result 
might come as a surprise: utility is often thought to be 
the socio-economic concept most similar to the physical 
concept of energy (or more precisely, the opposite of the 
energy), because agents seek to maximize their utility in 
social systems and physical particles minimize the energy 
in the zero temperature limit. Hence one might have in- 
tuitively expected the homogeneity of utility to be linked 
to a notion of temperature (the thermodynamic variable 
conjugated to energy), rather than to a chemical potential. 

1 K\\ equivalent formulation is to define the chemical potential 
\i as the logarithmic derivative of the partition function Z with 
respect to TV;, a definition that can be extended to some classes of 
noncquilibrium models 1141 . Note also that the standard definition of 
chemical potential for equilibrium systems differs by a conventional 
factor — 1/T from the one we use here [6]. 



Interestingly, Eq. (|16[) not only provides a link between 
two apparently unrelated concepts, but also yields a non- 
trivial prediction on the variations of utility across the 
system at non-zero temperature. As the chemical poten- 
tial remains uniform at any temperature, one sees from 
Eq. (fT5|) that the utility u q i = A, — Tds/dp q i becomes 
non-uniform if T > 0, and that the corrections to unifor- 
mity are given by the derivative of the local entropy. 

In a statistical physics language, Eq. (fT2|) corresponds 
to the canonical ensemble, where the number of interact- 
ing entities (agents or particles) is fixed. It is sometimes 
convenient to consider the so-called grand-canonical en- 
semble, where particles are exchanged with an external 
reservoir. In the context of agent-based models, the reser- 
voir corresponds to the external world. This means that 
we implicitly consider a very large set of sites ('the world') 
and focus only on a small subpart of it ('the system'), still 
containing a large number of agents. Since the 'world' 
has a fixed number of agents, it can be described by the 
stationary distribution Eq. ([S]) . Following standard statis- 
tical physics methods [5], the probability distribution of 
the considered subpart is given by 

V ow ({p ql }) = -L JJ e ff [/ 5 (^.-.^)-E i A iW ]/T ) (17) 
Zow q 

where Xi is the chemical potential of group i imposed by 
the external world. Finding the most probable densities 
p* qi is now straightforward since the densities on different 
sites are independent. Maximizing the argument of the 
exponential in Eq. ([T7|) . one recovers Eq. (fl6|) . 

In the following, we give two examples of models be- 
longing to the above generic class, in the context of urban 
economics. 

A simple urban economics model. The model 
presented here is a simple model of land use and transport 
interaction in urban economics |15j . In this model, a city 
is described as a grid composed of Q blocks. In each block, 
one or several agents (representing households) can live by 
paying a rent to the landowner. A central business district 
(CBD) is placed on the grid and all agents commute there 
for their work (monocentric city model) . A transport cost 
c per unit distance is associated to this commuting. The 
distance between a block q and the CBD is denoted by r q . 
The size of the city is fixed: a constant radius r / defines 
the urban fringe, out of which no agent livefl All agents 
have the same income Y, which is spent on transport, on 
housing and on a composite good z representing all other 
consumer goods. This gives a budget constraint for each 
agent 

Y = z + cr q + ap q (18) 

2 In standard urban economics models, land is used for agriculture 
outside the city, and the landowners then earn an agricultural rent 
1 1 5 1 - These landowners rent to the highest bidder, so that all prices 
must be greater than the agricultural rent. However, to simplify the 
presentation, we have dropped the agricultural rent parameter by 
introducing a fixed city size. 
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where a is the surface of housing, and p q is the rent per 
unit surface in block q. We first consider a simple model 
where all agents have the same surface of housing. Each 
block of the grid is composed of H cells of surface <jq. A 
configuration of the city is then given by the number of 
agents n q in each block q. We make the further simplifying 
assumption that the price p q of housing in a block q only 
depends on the density p q = n q /H of agents in this block, 
namely p q = p{p q ). Let us emphasize that this hypothe- 
sis is an important simplification with respect to standard 
urban economics models, in which the price emerges di- 
rectly from the competition for land between agents, and 
the density from their utility maximization with respect 
to the surface of housing |15j . In cases where an explicit 
expression is required, we will use a logarithmic form 

p(p q )=PoHl + p q ), (19) 

where po is a positive constant. 

The utility function has to be specified explicitly. It 
should be an increasing function U (z) of the quantity of 
composite good z each agent consumes, that we choose 
to be simply U(z) = z. This means that, in the limit 
T — > 0, each agent wants to maximize the share of his 
income which is left after transport and housing expenses. 
Using Eqs. ([T5jl and (fT§|) . the utility U becomes a function 
u q (p q ) of the local density, 

UqiPq) = Y - cr q - <j a p(p q ). (20) 

Urban economics distinguishes closed city models, 
where the total number N of agents is fixed, and open 
city models, where N fluctuates due to exchanges with 
the external world [IS] • We start by considering the closed 
city model. In the continuous limit where H and N — > oo 
with the average density ~p = N /(HQ) fixed, the station- 
ary probability distribution takes the form Eq. (|12[). with 
fq(Pq) g iv «i by 

fq{ P q)= U q (p)dp + Ts(p q ) (21) 

Jo 

and s(p q ) = -p q lnp g . 

The most probable density p* is then obtained as a func- 
tion of A from Eq. (|T5)) . namely 

Uq {p*)+T^-{p*) = \. (22) 
ap q 

In the limit T — > 0, often considered in socio-economic 
models, one finds p* — p*(r q ,\). with 

p*(r q ,\)=p- 1 ( Y - C ^- X y (23) 

where p" 1 is the reciprocal function of p. 

The parameter A is then determined from the density 
constraint p* = Q~p. Following standard literature 



[T5] . we focus here on the simplest situation of a one- 
dimensional city. Using the continuous approximation 

^E^fVM)^, (24) 

we compute the average density p(A), and then determine 
numerically the reciprocal function A(p). 

We now briefly turn to the open city model (similar to 
the above 'open world' case) where agents can also move 
to or from a large number of other cities. The stationary 
distribution is given by Eq. (fTT)) . which in the present open 
city model simplifies to 

r oc ({ Pqi }) = J-IJe ff ^-^/ r . (25) 

Zoc q 

Finding the most probable density is then an uncon- 
strained maximization problem. The relation df q /dp q = A 
yields the same equation as (|22p , resulting in the same den- 
sity profile ([23) in the limit T 0. For T > 0, the density 
can be obtained from a numerical resolution of Eq. (|22|l. 
We find that increasing the temperature T progressively 
blurs the zero temperature profile given by Eq. (|23[) . even- 
tually leading to a homogeneous density. The same effect 
has been observed in urban economics models [16] . As a 
consequence, the city is more spread, leading to a utility 
gain for agents near the city center, and to a loss for agents 
in the periphery. 

Urban model with two types of agents. In 

this second model, two income groups are distinguished. 
Rich agents (group 1) have an income Y\ and a surface of 
housing <7i, while poor agents (group 2) have an income 
Y 2 < Yi and a surface of housing a 2 < 0"i- Each block 
contains at most H agents, irrespective of their group. 
A configuration of the city is described by the densities 
Pqi = n q i/H and p q 2 = n q ijfL in each block q. The price 
P q i an agent of group i pays for housing in block q depends 
on his surface of housing and on the local density of poor 
and rich agents: 

Pql(Pql, P<fl) = V\P{aiPq\ + bip q2 ) ^ 
Pql{Pqr,Pq2] = 0- 2 p(a 2 p q l + b 2 p q 2) 

where a\, b\, a 2 and b 2 are given constants, and p a func- 
tion to be determined. The utility function of an agent of 
group i = 1, 2 in block q has the form 

u ql {p q i,p q2 ) =Yi- cr q - P qi {p q i,p q2 ). (27) 

The model is analytically solvable if Eq. (fTTj) is satisfied. 
For this condition to hold, one can choose a± = a 2 and 
b\ — b 2 . Then if o\b\ = u 2 a 2 , the function p can take any 
form, for instance the logarithmic form Eq. (|19[) used in 
the previous model, in which case we get (choosing b± = a 2 
and 0,2 = 01 ) 

Pqi(Pql i Pq2) = <*i PO ln(l + <J\P q \ +<J 2 p q2 )- (28) 
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Fig. 1: Density profile p qi as a function of r q for both groups of 
agents (rich, full line; poor, dashed line) for different temper- 
atures: T/Tq = 0.0018 (a), 0.018 (b), 0.089 (c) and 0.36 (d), 
with Tq = poo 2 . The dotted lines indicate the total density 
P* q i + P* q 2- Parameters: p = 1.4, c = 0.4, Yi = 452, Y 2 = 301, 
ci = 6, (72 = 4, Tf = 30, Pi = ]>2 = 0.13. 

The stationary distribution is given by Eq. (|12p. with 

fq(Pql,Pq2) = lq(Pql, Pql) + Ts{p qU p q2 ). (29) 



The expression of s(p q \,p q2 ) is given by Eq. (TI 
to = 2. Expressing l q (p q i, pqi) explicitly, we get 



with 



rPqi rp q 2 
lq{Pql, Pql) = \ U ql (p,0)dp+ U q2 (p q i , p) dp. 

Jo Jo 

(30) 

The validity of Eq. (jTDJ) , as well as the symmetry of 
Eq. (f3"U)) with respect to p q \ and p q2 , can be checked using 
Eq. (JTTJ). The equilibrium densities {p* q \,p* q2 ) are deter- 
mined from Eq. (|16[) . yielding a system of two non-linear 
equations, to be solved numerically. The results of this 
numerical resolution are shown on Fig. Q] One recovers at 
low temperature the standard separation, typical of north- 
american cities, between poor agents in the city center, 
and rich agents in the periphery [15| . The effect of a tem- 
perature increase is mainly to blur the zero temperature 
pattern, hence avoiding total segregation. 

Therefore, Eq. (fTB]) provides a direct prediction for the 
utility profile at arbitrary temperature T. It would be in- 
teresting to know whether this result remains valid beyond 
its a priori domain of validity, namely for models satisfying 
Eq. (fTTj) so that a function L can be defined. Considering 
again the above urban model with two types of agents, we 
keep the logarithmic form Eq. (|19[) for p, and choose as an 
example ai = a% = 1 and b\ = b 2 = 0. These values imply 
o\b\ ^ a 2 a 2 so that Eq. (fTT]) is not satisfied, ruling out 
the possibility to find a potential function L and to get a 
simple analytical solution of the model. 

Performing numerical simulations of this agent-based 
model with two income groups, in the case of a one- 
dimensional closed city, we first validate the simulation 



440 



400- 



360=* 



320 




280 



240 



200 



0.2 0.4 0.6 0.8 1 

r / r. 
q f 



160, 



1 i 1 i 1 i 1 i 1 



J i I i I i I lZ3 

0.2 0.4 0.6 0.8 1 

r / r, 

q f 



Fig. 2: Agent-based simulations of the model with two types 
of agents (left: i = 1; right: i = 2), when no potential func- 
tion L exists (one-dimensional closed city). The utility u q i is 
plotted as a function of r q for two temperatures (T = 100, *; 
T — 20, x). Adding to the utility the term Tds/dp qi yields 
the chemical potential Ai (T = 100, 0; T = 20, □), found to be 
constant throughout the system although the analytical solu- 
tion is not known. Full lines are horizontal lines indicating the 
spatially averaged value of the chemical potential. Parameters: 
Po = 70, c = 2, Yi = 480, Y 2 = 300, <n = 6.4, a 2 = 4, r f = 50, 
Ni = N 2 = 2000, H = 200, Q = 51. 



thanks to a comparison with the above solvable case. 
Turning to the non-solvable case, we test the validity of 
Eq. ([TO]), that is, whether the chemical potentials A.; = 
u q i + Tds/dp q i (i = 1,2) are uniform over the city for 
T > (when T —> 0, the utility should be uniform any- 
how). We indeed observe that for a non-zero temperature, 
the chemical potentials arc homogeneous even in this non- 
solvable model, while the utility is not (see Fig. [2J . We 
also note that although the number of agents is relatively 
small (H = 200), no systematic space-dependent correc- 
tion to the chemical potential is observed. 

We further use the obtained values of Ai and X 2 to per- 
form a numerical resolution of the system of non-linear 
equations (|16[) . which wc assume to be valid even in the 
absence of a potential L function. Interestingly, the results 
of the agent-based model and of the numerical resolution 
are found to be in very good agreement (see Fig. [3]), thus 
yielding a complementary test of the validity of the chem- 
ical potential approach. 

The validity of Eq. ([TB]) in this case can be understood 
as follows. In this paper, we focused on cases when the 
probability distribution has the factorized form Eq. (|12p. 
which is a consequence of the existence of a potential func- 
tion L. When no function L exists, the stationary distri- 
bution is no longer factorized, and we do not know its 
functional form. However, if the stationary distribution 
has only short range correlations, a chemical potential can 
still be introduced, in the same way as a chemical poten- 
tial can be defined in a physical system with short-range 
interactions [B]. 

This result is also consistent with a phenomenon known 
in the nonequilibrium statistical physics literature as 
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Fig. 3: Comparison of the density profiles obtained by agent- 
based simulations (i = 1, x; i = 2, o) and by numerical inte- 
gration (i — 1 , full line; i = 2, dashed line) using the values of 
Ai and A2 measured in the agent-based model (see text). Same 
parameters as in Fig. [2] with T = 20. 

"restoration of detailed balance" [TT1|T8] . Namely, start- 
ing from a microscopic stochastic dynamics which does 
not obey detailed balance, a coarse-graining procedure 
can lead to a detailed balance relation in terms of the 
effective, coarse-grained, degrees of freedom. This phe- 
nomenon can appear in physical systems in cases where no 
macroscopic flux is present (like fluxes of energy or parti- 
cles between boundary reservoirs), as macroscopic fluxes 
cannot be smeared out by the coarse-graining procedure. 
In the socio-economic models we consider in this paper, 
there is no obvious macroscopic flux, which suggest that 
detailed balance may be restored on large scales. 

Discussion. — In this letter, we have provided a 
clear relationship between the apparently unrelated no- 
tions of socio-economic utility and chemical potential. 
More specifically, we have shown that the uniformity of 
utility across the social system can be traced back to the 
conservation of the number of agents. This result not only 
provides a conceptually interesting link, but also yields 
non-trivial and testable predictions on the variations of 
utility among choices (e.g., sites, blocks) when T > 0. We 
also found numerical evidence that our result extends be- 
yond the class of models in which it was initially derived. 
It would thus be interesting to explore further its validity 
through numerical simulations of more realistic models. 

The idea of a non- uniform utility at equilibrium (Fig. [5J 
may be counter-intuitive for economists. Indeed, Nash 
equilibrium for homogeneous agents implies that all have 
the same utility, which seems not to be the case here when 
T > 0, since agents in the border of the city have a lower 
utility than those at the center. However, when noise is 
introduced in the decision process, a static equilibrium 
picture is no longer valid. Noise allows agents to explore 
the city, so that the time average value of utility is the 
same for all agents, leading to a macrostate described by 
Eq. (|16p through the ergodic hypothesis linking time and 
ensemble averages. The average utility of agents is then 
a decreasing function of T. Note that this picture of a 



"time-averaged agent" is close in spirit to the notion of 
"representative agent" advocated in discrete choice the- 
ory m|. It would be interesting to investigate further the 
relation between these two approaches. 

Another interpretation of our result is to consider the 
chemical potential Ai as an effective utility. We first note 
that the distribution V({p q i}) at T > can be obtained 
from the zero-temperature distribution by replacing l q by 
f q = l q + Ts [sec Eqs. (|T2j) and (|T3|) ]. in the same way 
as the macroscopic energy is replaced by the free en- 
ergy in a physical system at finite temperature. Then, 
changing l q into f q in Eq. (|TU| . we get an effective utility 
u qi = df q /dp q i. Hence the Nash equilibrium of an assem- 
bly of fictitious agents having this utility would precisely 
correspond to Eq. (|T6| , namely u q f = A, . 

# # # 

Fruitful discussions with S. Grauwin and F. Goffette- 
Nagot are hereby gratefully acknowledged. 
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